STRONG DIAMAGNETISM FOR THE BALL 
IN THREE DIMENSIONS 
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Abstract. In this paper we give a detailed asymptotic formula for the lowest 
eigenvalue of the magnetic Neumann Schrodinger operator in the ball in three 
dimensions with constant magnetic field, as the strength of the magnetic field 
tends to infinity. This asymptotic formula is used to prove that the eigenvalue 
is monotonically increasing for large values of the magnetic field. 



1. Introduction 
1.1. The operator and main results. Let fi be the unit ball 

fl = {x = (X!,X2,X 3 ) € R 3 | |x| < 1} 

and let B be a constant magnetic field of magnitude B > along the X3 axis, with 
a corresponding choice of magnetic vector potential A, 

B = (0,0,5), A=-{-x 2 , Xl ,0). 
We consider the magnetic Neumann Schrodinger operator 

H(B) — (— iV + A) 2 (1.1) 

with domain 

Dom(W(B)) = e W 2 > 2 (n) I N(x) ■ (-iV + A)*| an = 0}. (1.2) 

Here N(x) is the interior unit normal to dCl. This operator has compact resolvent 
and is semi-bounded from below, so it makes sense to enumerate its eigenvalues in 
an increasing order. For a self-adjoint operator TL that is semi-bounded from below 
we will use the notation \j,n to denote its jth eigenvalue. In particular, we will 
write 

^i,H(B) = inf Spec(H{B)) 
for the lowest eigenvalue of H(B). The first main theorem of this paper concerns 
the asymptotics of Ximb) as B — * 00. 

Theorem 1.1. There exist constants \j, j — 0, ... ,5, and Co; Cii C2; C3; $o> an d 
C such that with 

Cs(m,B)=m-|- CoVs - CiB 1/3 - (2B 1 / 6 . (1.3) 



and 



it holds that 



A B = inf CsK^-Cs (1.4) 



5 

Kn(B) = bJ2 X jB-^ 6 +S A 2 b +C + O{B- 1 ^), asB^oo. (1.5) 

]=0 



1991 Mathematics Subject Classification. 81Q10; 35PXX,82D55. 

Key words and phrases. Eigenvalue asymptotics, Large magnetic field, Unit ball, Ginzburg- 
Landau functional, Surface superconductivity. 

1 



2 



S0REN FOURNAIS AND MIKAEL PERSSON 



Remark 1.2. In the course of the proof we will obtain explicit expressions for the 
constants in the theorem above, especially it holds that 

A = Go, Ai = 0, and A 2 = 2- 2 ^0 S 1 /3 , (1.6) 

which agrees with the asymptotics of Xi^b) that is given in [14] for more general 
domains [/cl 3 (see Theorem II .4[) . The constants Go, A) an d #o are well-known 
universal constants which appear in the study of two model operators, see Appen- 
dixEl 

Before stating the next main theorem it is worth noticing that the constants 
Ao = Go and So in the theorem satisfy \ < Go < 1 an d < So < 1, so especially, it 
holds that 

0o-^o>O. (1.7) 

Theorem 1.3. Let So and Xq = Go be the constants from Theorem ll.il The 
directional derivatives 



Xl,H(B+e) — Xi.H(B) 
,H(B),± ~ £ ^ ± " 



K,n(B),± — 1™ 



exist and satisfy 



X 'i.-H(B),+ -Ku{B),-i forallB>0, (1.8) 
liminf X[ n(B) , > Go — -So > 0, and (1.9) 

limsupA' 1/H(B) _ < G + -S - (1.10) 

In particular, the function B t— > Ximb) is monotonically increasing for sufficiently 
large B. 

1.2. Motivation. 

1.2.1. Strong diamagnetism. Let Tlu(B) denote the magnetic Neumann operator in 
a bounded and smooth domain U C M n , n — 2,3. From the diamagnetic inequality 
(see [TB]) it follows that 

Xl,-Hu(0) < Xx t Uu{B) (1-11) 
for all B > 0. One might ask if the stronger monotonicity 

< Bi < B 2 X hHu{Bl) < X hHu{ B 2 ) (1-12) 

holds. In [4] counter-examples are given showing that (|1.12l) does not hold in general 
for all B\ and B 2 . The examples are given in R 2 , for constant magnetic field and 
the presence of a scalar potential, and for variable magnetic field without a scalar 
potential. 

Lately the question whether there exist a Bq such that Xi t Uu(B) is monotonically 
increasing for all B > Bo has been studied in detail. This is well-understood 
in two dimensions by now, with the final affirmative answer for regular domains 
in [71 [5j and for domains with corners in [3] . We discuss below the progress in three 
dimensions so far, which motivates our analysis for the ball. 

To continue, we introduce some conditions on the domain U. Let L c dU be the 
set of all points on the boundary where the magnetic field B is tangent to dU , i.e. 

T = {xedU | B-iV(x) = 0}. (1.13) 

Assumption 1. Let d denote the differential on dU . Then 

d(B • N(x)) ^ 0, for all xeT. 



STRONG DIAMAGNETISM FOR THE BALL IN THREE DIMENSIONS 



3 



If this assumption holds then T consists of a disjoint union of regular curves. 
We can orient them and denote by T(x) an oriented unit tangent vector at x G T. 
It is noted in [9] that this implies that the magnetic normal curvature k ny B{x) — 
K x [T(x) A N(x), p|r) is non-zero on T. Here K denotes the second fundamental 
form on dU . 

Assumption 2. The set of ■points in T where B is tangent to T is isolated. 

The following asymptotic formula of X x t Uu(B) was proved in [2] (the upper 
bound was given in [2"D]L 

Theorem 1.4. Let U C M 3 be a bounded and smooth domain that satisfies As- 
sumptions]]} and\^ Then there exist constants Go, 70 > and r\ > such that 

Khu(b) =QoB + %B 2 / 3 +0(B 2 ^), as 00. (1.14) 

The constant Oo is the same as in Theorem 11.11 and the constant 70 is given by 
70 = mf xer jo(x) where 

1 /3 

7o(x) = 2- 2 / 3 ,> 6 1 Q /3 \k n Mx)\ V3 (l + (So - l)\T(x) ■ it Q , 
and Sq and vq are fundamental constants given in Appendix [X] We note that 

70 = 2- 2 / 3 l >o<5o 1/3 (1-15) 

when U = Vt is the unit ball which makes our Theorem 11.11 compatible with Theo- 
rem [TT4J 

Using the expansion (11.14[) of the lowest eigenvalue \i,Hu(B)> the following mono- 
tonicity result was proved in [5] . 

Theorem 1.5. Let U C R 3 be a bounded and smooth domain that satisfies As- 
sumptions^ and\^ Let {ri,...,r n } be the collection of disjoint smooth curves 
making up T. Assume that for all j there exists x G Tj such that 70 (x) > 70. Then 
the function B > Xi _Uu(B) is strictly increasing for sufficiently large B . 

This shows that (|1.12|) holds for large values of B\ and B^. Even though the 
Assumptions[T]and[2]are fulfilled for the ball, the assumption on 70 in Theorem ll.5l is 
not. Indeed, for the unit ball O, the set T consists of the equator {x G dfl \ x% = 0} 
and the function 70(2;) is constant 

7o(z) = 2- 2 / 3 i>o<5o /3 , ier. 

1.2.2. Superconductivity. We consider superconductivity in the Ginzburg-Landau 
model. For a superconducting material of shape U subject to an external magnetic 
field Ka[3, with (3 — (0,0, 1), the Ginzburg-Landau energy functional is as follows, 

a) = / + naa)^ 2 - K 2 \V\ 2 + ^-|*| 4 cLr 

Ju z 

+ {no-) 2 j | curia -/?| 2 da;. (1.16) 



Here k > is a material dependent parameter called the Ginzburg-Landau param- 
eter. For given k, the parameter a measures the strength of the external magnetic 
field. The function ^ G i/ 1 (?7) is in this context called an order parameter and 
1 measures the local superconducting properties (density of Cooper pairs) of the 
material. Finally a G ff^QR 3 ) is the induced magnetic vector potential. In order 
to get a well-defined minimization problem, one uses gauge invariance to restrict to 
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vector potentials satisfying diva = 0, and impose the finite energy condition (see 
Q3J] for details), 



/ | curia -fl'd*<oo. 



A state a) where \t = and curia = j3 is called trivial. 

The analysis of magnetic ground state eigenvalues described above is relevant 
in superconductivity for the understanding of the loss of superconductivity in the 
presence of strong magnetic fields. The value of the magnetic field strength at 
which the material loses its superconducting properties, i.e., the minimizers of the 
Ginzburg-Landau functional have \& = 0, is called the third critical field. The 
calculation of this critical field has a long history, see [TTl[Tl [TTl[Tgl[T!)l lTllgll5]. 

In [9l [5] it was proved that, for k sufficiently large, the following sets are equal 

jV"(k) :— {a > | £ Ki( j has a non-trivial minimizer} 

= {cr > | £ Kt „ has a non-trivial stationary point} (1-17) 

= {a > | Xi t nu( K a) < k 2 }- 

In the case where U is the unit ball f2, we can use the monotonicity result of 
Theorem 11.31 to conclude that, for n sufficiently large, 

{o0 | A 1)Wn(KCT) < « 2 } = (0,HcM)- (1-18) 
Here a — Hc 3 (k) is the unique solution to the equation 

Hereby, we get a complete determination of the third critical field, Hc 3 (k), for large 
values of n. Upon inserting the asymptotics (|1.5[) . one gets a six-term asymptotic 
expansion of Hc 3 (n) for the ball, 

e e^ /3 &l' 2 \3e 7 /3 e 4 > 3 J V6eJ 3/6 e 7 /6 , 

A| A470 _ j|_ _ a Ihc 3 («) 

2eg eg 3eg e 

Thus, SS(k) is an interval (for large k) both in the case where U satisfies the assump- 
tions of Theorem 11.51 and in the case where U is a ball. It remains an interesting 
open question to prove this result in general, i.e., to prove strong diamagnetism for 
general (smooth) domains in R 3 . 

1.3. Organization of the paper. The main part of this paper is devoted to the 
proof of Theorem ll.il The proof is divided into several parts: 

We denote by 7i m (B) the operator TC(B) restricted to angular momentum m € Z. 
In Section [2] we show a lower bound for 7i m (B). First we show in Lemmas 12.71 
and 12.81 that (in terms of m and B) either the first eigenvalues of TC m (B) are too 
large to be compatible with (jl.5p . or we can reduce 7i m (B) to an effective operator 
Q m (B), satisfying, as B — * 00, 

hn m (B) = B\ 3 . Qm[B) + 0(B^ 2 ), for 3 = 1, 2. (1.19) 

For the values of m and B such that (|1.19p holds, we prove a lower bound for 
Q m (B) in Proposition EH 

In Section[3]we use the lower bound from Scction[2]to give localization properties 
of eigenfunctions of Q m {B). These are used in Section 2] to obtain a spectral gap 
formula for Q m (B), which together with (|1.19[) implies a spectral gap formula for 
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7~Lm(B), showing that for some 7 > it holds that (still under some restrictions on 
m and B) 

\n m {B) >e o B + (7o + 7 )5 2 / 3 + 0(B 7 / 12 ), asS^oo. (1.20) 

In Section [5J Theorem I5.1f iii) we use the Grusin method to, for certain m and 
B, calculate a trial state that together with the spectral gap formula gives upper 
and lower bounds on \\^u m (B) which are compatible with (|1.5[) . We also give two 
alternative trial states, in Theorem 15. lf i) and (ii), for values of m that are further 
away from the optimal choice. 

In Section[5]we show that X 1 t n m (B) is larger for the values of m and B not treated 
in Theorem I5.1f iii). Depending on m and B we use different methods to achieve 
this. For m and B which are far from the optimal region we use (a refined version 
of) the lower bound from Section [21 For m and B that are closer to the optimal 
region we use the trial state from Theorem I5.1f i) and (ii) which by the spectral 
gap formula (|1.20p must be \\,H m (B)i but which is strictly greater than the ones 
obtained in Theorem l5. If iii) . Finally, in Section[7]we minimize the eigenvalue found 
in Theorem 15 . If iii) ■ which proves Theorem ll.il 

The proof of Theorem ll.3l is a consequence of Theorem ll.il using a perturbation 
argument from [5]. The details are given in Section [5J 

We start by introducing the new coordinates and some quadratic forms and 
operators. 

1.4. New coordinates, auxiliary operators and quadratic forms. We con- 
clude this first section by introducing the coordinates we will work in and the dif- 
ferent quadratic forms we will work with. First, we switch to spherical coordinates 
(£1,2:2,2:3) >-> (r,(fi,0), 

x\ — r cost/? sin 

x 2 = r sin sin 0, < r < 1, < 9 < n, < tp < 2n. 
X3 = r cos 9 1 

We decompose the Hilbert space as 

L 2 {VL) = L 2 ((0,1) x (0,7r),r 2 sin(9drd6l) <g> L 2 {S 1 ,dip) 

-im<p 



- ^ 2 ((0, 1) x (0,7r),r 2 sin6»drd6») 



m— — 00 

2/ 



/2tt 

that is, for a function ^ E L 2 (fi), we write 



'2tt 



where ip m G i 2 ((0, 1) x (0, 7r), r 2 sin#dr d9) . Next, we write the operator TL(B) 
corresponding to this decomposition as 

00 

H{B)= H m {B)®l 



m— — 00 



where H m (B) is the self-adjoint operator acting in L 2 ((0, 1) x (0, ir), r 2 sinfldr d^), 
given by 

n m (B) = -o 2 - 2 -d r -\d 2 -^- de /B, '^° 



r ' r 2 tan 9 V 2 r sin ( 

with Neumann boundary condition at r = 1. In the continuation we skip the 
subscript m on ij) m and write just ip. Inspired by [14] we introduce the new scaled 
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coordinates (r, p) as 

( T = VB(l-r), 

< " (1-21) 

lp=^-7r/2), 

with corresponding new domain 

0(B) = {(r,p) | < r < VB, -\VB < p < \Vb}- 



In fact, for a point x in f2, r = \AB dist(x, <9il) is equal to the (scaled) distance to 
the boundary and for a point x G dfl we have p = \fB distant, T), the (scaled) 
distance along the boundary to the equator. 

The quadratic form corresponding to j^Ttm (the prefactor 1/B is just for conve- 
nience) transforms into the quadratic form 



n(B) 



2 



1 ( B(l-B-^T)cos(B-^p) m_ _ _ 

B\ 2 (l-B-VVJcos^-VSp)^ ? 

x (1 - £T ^V) 2 cos(B- 1 /'V)S- 5/6 dr dp 
in the Hilbert space 

L 2 (fi(B), (1 - B-^V) 2 cos(B- 1 / 3 p)B- 5 / 6 dr dp 

We apply the unitary transform Hip = B~ 5 / 12 ip to get rid of the factor £>~ 5 / 6 in 
the measure and work in the Hilbert space 

L 2 (n(B), (1 - B- 1/2 t) 2 cos(B- 1/3 p) dr dp) 

instead. We will, by abuse of notation, continue to write ip instead of Utp. We 
denote by Q m {B) the operator corresponding to the quadratic form q m . 

Next we want to define a quadratic form q m by the same integral expression as 
for q m but in the Hilbert space L 2 (R+ , dr dp) , where 



3 2 



{(t, p) I < t < oo, — oo < p < oo}. 



However, since neither 1 — B~ x I 2 t nor cos(B _1 / 3 p) are strictly positive in K 2 ^, we 
make a technical modification to be able to talk about the corresponding operator. 
We define smooth functions I : R+ — > R+ and Cos : R — > M that satisfy 

„/ x I £ < 4 5 , ^ . . I cos(x), \x\ < ?, , „„. 

l(x) = \ 1 > and Cos(x)= , U ' " (L23) 

^ 2 ' — 2 ' ^ 2 ' I 2 "! — 3 ' 

and such that I is monotonically increasing in the interval (5,5) and Cos is even 
and monotonically decreasing in the interval (f,f)- The quadratic form q m is 
defined by 



MA = [ 



2 B" 1/3 2 
l^f + " , -— o 1^1 



(l-^fl-VV))' 
1 / J B(l-£(S- 1 /2 r ))Cos(B- 1 /3p) 



2 



+ B V 2 (1 - ^(-B -1 / 2 t)) Cos(S- 1 /3 p ) J 1^1 

x (l-^(B- 1 / 2 r)) 2 Cos(B- 1 / 3 p)dTd /9 . (1.24) 
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We denote by Q m (B) the self-adjoint operator that corresponds to q m . It is an 
operator in L 2 (r 2 + , (1 - ^fl-VV)) 2 CosCB-VSp)) with Neumann condition at 
t = 0. An integration by parts, show that it acts as 



, 2£'(B~ 1 / 2 t)B~ 1 /' 2 



Q m (B) = -8< r+ { ;_ i(B . 1/2T)) 9 T 



B- 1 / 3 (a2 Cos' \B- 1 ^p)B- 1 ' 3 



J 2 uos^g r 



(l-/(B-Va T )) 2 V " Cos(B-V3 p ) 
1 - 1{B- x I 2 t)) CosiB-^p) 



B\ 2 (1 - ^(B-Va T )) Cos(fl-V3 p ) 

Finally, we also define the quadratic form q m in L 2 (E^_, drdp) by 



(1.25) 



<b»bl>]= [ \d T ^\ 2 +(r + ^={m-B/2) + B- 1 l^\ \i»\ 2 + B^ 3 \d p ^\ 2 dr dp 
Jr 2 + \ VB 2 J 

(1-26) 

with corresponding self-adjoint operator Q m (B) with Neumann boundary condition 
for r = 0. 

2. A ROUGH LOWER BOUND 

In this section we recall the localization formulas of the lowest eigenfunction of 
H(B) obtained in [Uj and written out in detail in [5], use them to reduce the study 
of Tim(B) to the study of Q m (B), and give a rough lower bound of its quadratic 
form q m in Proposition ^. 91 

Theorem 2.1. Suppose that U C M 3 satisfies the Assumptions]]] and\^ and that 
<f satisfies H{B)^ = A* with X(B) < Q B + luB 2 ^ 3 . Then there exist positive 
constants a\, a%, do, C and Bq such that 

f e 2a x B^ dist(x,9Q)^|2 + £-l|(_ iV + A)4-| 2 ) dx < C||*|| 2 

Ju 

and 

f e 2a 2 BV2 dista[/(Xir) 3/ 2 (|^|2 +jB -l| ( _. v + A) ^|2- ) dj; 

Jidist(x,dU)<d n \ 



for all B > B . 



<C7e CBl/8 |!vI/|| 2 



From now on we assume that U — f2 is the unit ball. Then the set T, intro- 
duced in (|1.13p . consists of the equator and we can extend the distance function 
distant, T) to all of SI (except the origin) as distant, T) = distao(i,r) where 
x = t^t S dfl. By the exponential decay away from the boundary, the second in- 
equality in Theorem 12.11 is then valid with the integral on the left-hand side being 
over all of fl, with possible changes of the constants. We will use the following 
corollary. 

Corollary 2.2. Suppose that * satisfies H(B)^> = A* with \(B) < Q Q B + ujB 2 / 3 . 
Then for all n € N there exist positive constants C n and B n such that 

f dist(x,ari)"(|*| 2 + J B- 1 |(-iV + A)*| 2 ) d.x < C7„B-™/ 2 ||*|| 2 (2.1) 
Jn 

and 

f dist an (a;,r)"(|*| 2 + B- 1 |(-zV + A)*| 2 ) dx<C7„B^ l/4 ||*ir (2.2) 
Jn 
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for all B > B n . 

Remark 2.3. The order in (|2.2p is not optimal. The calculations below indicate 
that the same estimate is true with B~ n l 3 instead of i?~™/ 4 in the right-hand-side. 

Let < e < 1/12 be given. We introduce a smooth cut-off function < xb < 1 
such that 

jl, if dist(x,dfi) < B- 1 / 2 ^ and dist 9n (a;,r) < £? -1 / 4+E , 
XB ~ [0, if dist^,^) > 2 J B- 1 /2+e or dist^^^r) > 2S- 1 / 4+e , ' " " ' 

and such that |Vxs| < CB 1 ! 2 ^ for some C > 0. 

Lemma 2.4. Suppose that * satisfies U{B)^ = AW iwi/i A(B) < 6 S + wB 2 ' 3 . 
For any N > it holds that 

[ |(-iV + A)*| 2 da; = f |(-iV + A)(xs*)| 2 da; + C'(S~ Ar )||*|| 2 , as B -» oo. 

Proof. This follows by commuting (— iV — A) and \b and using Corollary 12. 21 □ 

We remind the reader of the quadratic forms q m and q m , introduced in (|1.22[) 
and (|1.24p . with corresponding self-adjoint operators Q m (B) and Q m (B). First we 
note that 



inf Spcc(H{B)) = inf inf Spec(H m (B)) = B inf inf Spec(Q m (S)) . (2.4) 

We use Lemma [2.41 to reduce the study of H(B) to the study of the quadratic 
form q m in the half-space R 2 ^. We will denote by ip B the function 

1pB = XBlp, 

where \b is the cut-off function from (|2.3|) . Notice that ^, and by consequence ip B , 
depends on m, but we do not include this in the notation. 

Lemma 2.5. Assume that ipi and ip2 satisfies Q m (B)tpj — Xj(B)ipj with Xj{B) < 
OqB + loB 2 / 3, . For any number N > there exist constants Bn and Cn (indepen- 
dent of m) such that if B > Bn then 

|q m (^i,^2) - q m (XBtpi,XBip2)\ < C N B- N \\ipx\\ ■ \\ip 2 \\- 

Proof. This follows from Lemma 12.41 transforming to the coordinates (r, p) . □ 

We recall that < e < 1/12 and note that 

suppV-s C {(t, P ) I < t < 2B £ , ~2B 1 / 12+£ <p< 2B^ 12+£ }, (2.5) 

which implies that B~ 1/2 t < 2B £ ~ 1/2 and \B~ 1/3 p\ < 2B £ ~ 1/4 on the support of 

IpB- 

Lemma 2.6. Assume that ip satisfies Q m (B)ip = X(B)ip with X(B) < QqB + 
ljB 2 / 3 . For any number n > there exist constants B n and C n ( independent of m) 
such that if B > B n then 

[ r n (\^ B \ 2 + \d T ^ B \ 2 + B^ 3 \d p ^ B \ 2 )dTdp<C n U B \\ 2 (2.6) 

and 

I \p\ n {\M 2 + \drM 2 + B- 1/3 \d p1 p B \ 2 )drdp<C n B n / 12 U B \\ 2 . (2.7) 

Jm 2 + 

Proof. This follows directly from Corollary 12.21 by transforming into the new 
coordinates (r, p). □ 
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We will use the estimates in Lemma 12.61 to reduce the values of the angular 
momentum m that we must consider, which in the end will enable us to study the 
effective quadratic form q m instead of q m . Let Co = -^g( m ~ B/2). 

Lemma 2.7. Assume that if> satisfies Q m (B)ip = X(B)ip with X(B) < QqB + 
u>B 2 / 3 . There exist positive constants D and Bq such that i/|Co| > D and B > Bo 
then 

* m \M > ^CoIIV'bII 2 . (2.8) 

Proof. We expand the potential in q m and collect the terms in front of the different 
degrees of Co- By (TmH and we see that if B > max((8/ 7 r) 4 /( 1 - 4e ), 6P/(i-2e))) 
then it holds that 1{B~ x / 2 t) = B~ x I 2 t and Cos(B' 1 / 3 p) = cos(B~ 1 / 3 p) on the 
support of ijiBi an d so we can write the potential in q m as 

'v A B(l-B- 1/2 r) 2 cos 2 (B- 1 /3 /5) /B\l 2 1 



Co 



cos(B- 1 /3 p ) 



2 2 
Using that < cos(i?~ 1 / 3 p) < 1 on the support of ipB and the general inequality 

3 T 2 _ 1 
4 X 3 



(x — y) 2 > \x 2 — \y 2 , valid for real x and y, we find that qm^s] is bounded from 



below by 

1 /VB(l-S- 1 /2 T )2 cos 2 (B -l/3 j0) I 



4 S0 3 V 2 2 J cos(B-V 3 p) 

Using (|2.6p and (|2.7p we get the existence of a constant C such that 



\iPb\ 2 dr dp. 



n/B(1-B-VV) 2 cos 2 (B-V3 p ) y/sy 1 



cos(S- 1 /3 /9 ) 



This clearly implies ff278[) . □ 

By Lemma 12.71 above we need only to consider bounded Co ■ This enables to 
study the quadratic form q m instead of q m . We assume that |Co| < Do for some 
Do > and also let D = D + |Co| so that the inequality | Co — Co I < A) holds. 
Here Co = Co is the constant from Lemma IA.21 

Lemma 2.8. Assume that ipi and ip2 satisfies Q m (B)'ijjj — \{B)^jj with Xj(B) < 
QqB + ujB 2 / 3, . If |Co| < Do for some constant Dq > then there exist constants 
C > and Bo > ( independent of m and B ) such that 

\q m (XBipuXBi>2) - q m (XB^i,XBih)\ < C5- 1/2 ||V>i|| • U2W 
forB> B . 

Proof. This follows by expanding the terms in q m and estimating using (|2.6|) 
and ([277]) . □ 

Next, we introduce 

Ci = (Co - Co)S 1/6 = (^=(m - 5/2) ~ Co)s 1/6 , (2.9) 
and note that if | Co I < -Do then 

lCi|<A)B 1/6 - (2.10) 
With this notation, the form q m reads 

* m[lP] = L l^| 2+ ( T + ^ + ^ 1/6 (^ 1 + ^)) H'\ 2 + B- 1/3 \d p ip\ 2 drdp. 
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Given constants C\ , C 2 and M we define the function Wb as 

w (\ /So + CxB- 1 / 3 , \p\<M, 
Wb{p) = < „ u4n _ 1/3 !sl , (2- 11 ) 
[<d + C 2 M 4 B |p| > M. 

Proposition 2.9. Assume that |£i| < DqB 1 ^ 6 . There exist positive constants C\, 
C 2 , M and B such that for all rn G Z 



q™M > / W B (p)|V| 2 drdp 

JB2 



/or a// B > B and -0 G Dom(q m ) i/ M > M. 

Before proving this proposition we emphasize that these estimates do not only 
hold for ground states, but for all functions ip in the domain of q m . 



Proof. We will prove a slightly stronger statement than the one in Propositionf 
Let D > be a large number, to be specified below. In the following lemmas we 
consider the three cases: 

(1) Ci > D as B > Bi for some B x > (Lemma |2~TU1) . 

(2) Ci < —D as B > B 2 for some B 2 > (Lemma 12. 111) . 

(3) |Ci| < D for all B (Lemma EH]). 

In fact, for the cases (1) and (2) we will prove stronger estimates. □ 

Lemma 2.10. Assume that £1 > D and that B > D 6 for some D > 0. Then there 
exists a positive constant C , independent of D, such that 

q m [V] > (©0 + CD 2 B- 1 ^) U\\ 2 (2.12) 

for all tp G Dom(q m ). 

Proof. By the assumption on £1 we have 

3- 



for all Using the de Gennes model operator Q from Appendix IA.2I we get a 

positive constant C such that 

Ai,g (Vo + £T 1/6 (Ci + y) ) > A i,s (Co + DB- 1 / 6 ) 

>e a + CD 2 B- 1/3 , B>D 6 , 
from which ([2~T2]) follows. □ 

Lemma 2.11. Assume that |Ci| < DqB 1 ^ . Then there exist positive constants C 
and D such that if B > Dq S , D > D and Ci < —D then 

2 



q m M > (Qo + cd^b-^uw 



for all ip £ Dom(q TO ). 



Proof. We assume that Ci < —D, for some constant D > 0. Along the proof we 
will get some constraints on D that finally will determine D. 

We first assume that D > 1. Let < Xi,b{p) < 1 be a smooth cut-off function 
that satisfies: 

(A) xxAp) = 1 if (1 - ^ICi|- 1/4 ) v^l < \P\ < (! + ^ICi|- 1/4 )v5i]- 



(B) XliB (p) = if |p| < (1 - |Ci|- 1/4 )V^i or \p\ > (1 + |Ci|- 1/4 )V2|Ci|- 

(C) Ixi b(p)| < 'llCiT 174 for some constant ?i > 0. 
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(D) The function X2.b(p) 
some constant I2 > 0. 
The IMS formula gives 



I-xIb(p) satisfies \x' 2 , B {p)\ < ^2 |Cx |~ 1/4 for 



2 f 2 

q™W = q m [xi,Bip} + qmlx2,Bip}- b- 1 / 3 ^ / W. B { P )\ 2 • \^\ 2 drdp 



The error term above is bounded as 
2 



^ 1/3 E / |x^(p)r-M 2 drdp<(^ + ? 2 2 )icir l/2 s- l/3 nvii 2 

j=l 

From (B) we see that the support of xi,b is included in the set 

{p I Ci + ^>(-2|Ci| 3/4 + ICi| 1/2 )orCi + ^<(2|Ci| 3/4 + lCir /2 )}, 
so especially if |d| > 16 it holds that 



(2.13) 



5 + £ 



<fic.r 



on the support of Xi,b- By the estimate |Ci| < DqB 1 /^ we have 

^>2 



6 + £ 



- 2 



on the support of xi,b- If B is sufficiently large (B > £>q 8 ) we can Taylor expand 
Ai.g to find a positive constant C such that 

Kg (Co + B- 1/6 (Ci + ^)) > e + c(Ct + y) V 1/3 
for all p on the support of xi,b- We insert this into q m , to get 

2 2 

e |xi,s^| 2 + S- 1/3 (|9 P (xi, s V)| 2 + C(Ci + y) IXi.sVf)] drdp 



> 



> 



(e + C 1 / 3 2- 2 / 3 A 1 , A4 (^321/3^)5-1/3) || Xl sV) || 



where we have used the notation of the Montgomery model in Appendix IA. 31 By 
Lemma [A. 71 it follows that if |£i| is sufficiently large then 



A 1 ^(C 1 / 3 2 1 / 3 Ci) > 2C 1 / 6 2 1 / 6 |Ci| 1/2 , 



and so 



ImbCiM] > (©o + 2 1 / 2 C 1 / 2 |Ci| 1/2 i?" 1/3 )||xi^l| 2 

> (Go + 2 1 /2 C 1 /2^1/2 B -l/3 ) \\ XltBl pf. 

In the same way we show that on the support of X2,b it holds that 



(2.14) 



Ci + 



> ^ICil 3/4 



if |Ci| > 1/16. This implies that there exists a constant C > such that 



1/6 



\X2,B^\\ 

>(eo + qCil 3/2 ^ 1/3 )llx2,B^ir 
>(e a + CD^B-^)\\x2^A\ 2 - 



(2.15) 
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The proof is completed by combining the equations (|2.13[) . (|2.14[) and (12 . 1 5[) for a 
sufficiently large D. □ 

The case left to study is when there exists a constant D such that |£i| < D. 

Lemma 2.12. Assume that |£i| < D. There exist positive constants C\, C%, M 
and Bo such that with Wb from (|2.11|) 



q m [ip] > I W B (pM 2 dTdp 



for all B > B and ip € Dom(q m ) if M > M. 

Proof. Fix M > 0, to be specified below. Let us introduce a smooth cut-off 
function Xi,M that satisfies the following properties 

(i) 0< X i,m(p) < 1 

(ii) xiMp) = 1 if \P\ < M 

(iii) xiMp) = if \p\ > 2M 

(iv) There exists a constant l\ > such that \\\ m(p)\ — h/M for all p. 

(v) The function X2,m{p) = v/l~ Xi,m(p) satisfies |x 2 ,mG°)I ^ h/M for some 
constant I2 > 0. 

By the IMS formula, it holds that 

qtnbP] = qm[Xl,MV>] + <\m[X2,M^] ~ \ 2 I^,m(p)| ' \tp\ 2 dp. 

3=1 jR l 

The localization error is bounded by 

B ~ V3 E I ^m(p)\ 2 -W^dp<{ll + ll)^B-^\\^ (2.16) 

i=i Jr2 + 

so by choosing M large, we can make this error small. We locally introduce the 
notation tpj = Xj-M 1 ^, for j = 1,2. On the support of ipi we can use the Taylor 
expansion of Xi.g to get the existence of a constant C > such that 

Ai, e ( Co + s- 1/6 f Ci + ) > e + cs-va ft + 9 



for £? large enough. This gives 
qm[V"i] > 



^ 6o|^i| 2 + B-^[\d P M 2 + c(Ci + y) Vl 2 ) drdp. 

Next, we use the modified Montgomery operator in Appendix IA. 31 to get 

q B [W> / eo|^i| 2 + J B- 1/3 C 1 /3 2 - 2 /3A 1 , A4 (C 1 /32 1 /3^ 1 )| V;i |2 dTd/9 . 



Since Ai,^4 has minimum uq we get a constant C\ = ^C 1 / 3 2 2 / 3 Ai such that 

q«,W > (e + 2C 1 B- 1 / 3 )j|^i|| 2 
On the support of ^2 we have 
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where M is chosen large enough, so that the last inequality holds for M > M. 
There exists a constant C 2 > 0, independent of M, such that 

Ai, e (Co + 5- 1/6 (Ci + y ) ) > Ai,s (Co ± B^Hl 2 /10) 

> 9 + 2C*2M 4 B- 1/3 

and so we get 

c\ m [H > (O + 2C 2 M 4 B- 1 / 3 )\\i> 2 \\ 2 . 

This finishes the proof if we choose M so large that the localization error (|2.16p is 
dominated by d x / 3 1 1 -0 X 1 1 2 + C 2 M i B- 1 / z \\^ 2 \\ 2 . □ 

3. An improved localization formula 

Proposition 3.1. Let lu > and a > 0. Then there exist positive constants Bq 
and Cq such that if B > Bq and ip satisfy Q m (B)4> = AV> with A < 0o + wB -1 / 3 
then 

e 2 °H(|V;| 2 + |9^| 2 + fl- 1 /3|^|2) dTdp < Co || V; ||2_ 

"■+ 

Proof. Let xi.m and X2,m be the cut-off functions from the proof of Lemma \2. 121 
where M is going to be specified below. Also, for e > 0, let 

This function w £ is bounded and continuous on K and differentiable everywhere ex- 
cept at 0. Moreover |?4(p)| < 1 for all p ^ 0. In particular the function X2,M^e aVs ^ 
belongs to the domain of q TO . We use integration by parts (the IMS formula) to get 

X\\X2,M^ aVelp) \\ 2 = ^m[x2,M^ aMp) ] ~ I \d P (x2,Me av ^)^\ 2 dTdp. (3.3) 

Next, we choose M such that both (C 2 M 4 -u- 2a 2 ) > 1 and M > M hold, where 

M is the constant in Proposition 12.91 and we choose B greater than the constant 
Bq in Proposition ^. 91 Using the assumption on A and the lower bound on q m from 
Lemma \2. 121 we get 

(C 2 M 4 -u- 2a 2 v' e {p) 2 ) \ X 2,M^ avAp) T dr dp 

<2 / \ X '2,M(p)\ 2 e 2aMp) \i>\ 2 dTdp. (3.4) 

JUL* 

The function x' 2 m(p) * s supported in the set {p £ K | M < p < 2M}, where also 
the inequality e al,e< - p - 1 < e 4aM holds. Inserting this and the choice of M in (|3.4|) we 
get 

/ |x2AfV'e a, ' e(rt | 2 drdp< 2e 4aM / l^fdrdp. (3.5) 

Since the right-hand side is independent of e we can let e tend to zero and use 
monotone convergence to get 

/ X 2 2M z 2a]p] m 2 dTdp<2e iaM f \^\ 2 dTdp. (3.6) 

Jr 2 + Jr 2 + 

Since xi.m is supported in {p 6 K | |p| < 2M} we have 

/ X?,Me 2a|p| |^| 2 drdp<e 4oM / M 2 drdp. (3.7) 

JR 2 , 
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Combining (|3 . 6|) and (|3 . 7[) gives the L 2 -bound in (f3~T]) . Next we turn to the terms 
involving derivatives. Using the triangle inequality, we have for the p-derivative 

2 
+ 

<2 [ \d p (e av ^ X 2.Mip)\ 2 dTdp + 2 [ \d p (e a ^ X2M )4>\ 2 dTdp 

<2B 1 / 3 q m [e av ^ X 2,M^] +2 / |d p (e<"^ X 2.M)^| 2 drdp. 

The corresponding inequality for the r-derivative is, since v e and x?,m are inde- 
pendent of t, 

/ e 2av ^xl M \d T ^\ 2 drdp < q m [e av ^ X2 , M ^]. 
Jr 2 + 

Combining these two inequalities with (|3.3p gives 
3 2avAp hl M (\d^\ 2 + B- 1 / 3 ^ 2 ) drdp 

<2B~ 1 / 3 [ |9 P (e Q ^ ( ' 5) X2.M)V'| 2 dTdp + 3q m [ e a ^^X2.M^ 
Jm 2 + 

<(2B-^ 3 + 3) [ |9 P (e a ^ ( ' ,) X2,M)^| 2 dTdp + 3A||x2,MV'e Q ^ (p) || 2 . 



e 

2 



Moreover, 

\d P {e av ^ X2M )\ = K(p)+X 2 ,Af(p)|e aMp) < (« + h/M)e av 'V>, 
so for B > 1 we can use (|3.5p to get 



e 2av ^xl M (\d^\ 2 + B- 1 / 3 ^ 2 ) drdp 

< (5(o + / a /M) + 3(9o + w)) ||x 2 ,AfV>e a,ap) || 2 

< 2[5(a + ; 2 /M) + 3(e + w)]e 4aM / \^\ 2 dr dp. 

Jm.1 

By monotone convergence we have 

e 2alpl X 2 2M (l^| 2 + B- 1 ' 3 \d p ^\ 2 ) dr dp 

< 2[5(a + / 2 /M) + 3(e + w)]e 4aM / |^| 2 drdp. (3.8) 

The same estimate with Xi,M m place of X2,m is easier since we do not have to use 
v £ and the functions involved have compact support. The result is 

e 2a ^xl M (\d^\ 2 + B- 1/3 \d p ^\ 2 ) drdp 



< 



2[5(a + /i/M) + 3(6 + w)]e 4aM / (i/fdrdp. (3.9) 



Finally, a combination of the equations (|3.6|) , (|3.7p , (|3.8p and (|3.9p implies (|3.ip . □ 
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Corollary 3.2. For all lo > and neN there exist positive constants B n and C n 
such that if B > B n and ip satisfies Q m ip — with A < 0o + wB -1 / 3 then 

\ P \ n {\M 2 + \driJA 2 + £r 1/3 |d p V| 2 ) drdp < C n U\\ 2 . 

Proof. This follows from Proposition ^. 11 by noting that all the terms in the Taylor 
expansion of e 2a ' p ' are positive, and thus for all non-negative integers n it holds that 
e 2a \P\ > (2a\p\) n /n\. □ 

4. Improved lower bounds and a spectral gap 

Proposition 4.1. Let Ci be as in ([21?]) . i.e., (l = (^( m - B /2) - Co)B 1/e and 
X 1 ^ be the lowest eigenvalue of the Montgomery operator, see Avvendix \A.3[ 

(A) For all G\ > there exist constants Bq and Cq (independent of m) such 
that if \Ci\ < Ci and B > Bq then 

q m [V] > {Qo + \m{Ci)B-^-C B-^)\\^ (4.1) 

for alli/j S Dom(Q m (B)). 

(B) For all C\ > i/iere ea;is£ positive constants j, Bq and Cq (independent of 
m) such that if |Ci| < C\ it holds that 

\n m (B) > ®oB + (70 + 7 )S 2/3 - C Q B 7/12 (4.2) 

if B > Bq. In particular, if £i is bounded there exists a positive constant 
B\ ( independent of m ), for B > B\ then the set 

Spec(H m (B)) n (-oo, Q B + (% + i/2)B 2 ' 3 ) 

is either empty or consists of the lowest eigenvalue ofH. m (B). 

Proof. We recall that 



J r2 M 2 + ( T + Co + S^ 1/6 (Ci + y)) H'\ 2 + B- 1/3 \d p iP\ 2 dr dp. 



We start with the proof of (A). Fix < <r < 1/12. Let us introduce a smooth cut-off 
function < Xi,b{p) < 1 that satisfies the following properties 

(i) xiAp) = 1 if \P\ < 

(ii) X i,b(p) = if \p\ > 25 s 

(iii) There exists a constant l\ > such that \x'i b(p)\ — hB~* for all p. 



(iv) The function X2,b{p) = y 1 — x\ b(p) satisfies \x' 2 b(p)\ — hB~ q for some 
constant l 2 > 0. 

We denote by ipj = Xj,Bi\>- Then clearly both ipi and ip2 belong to the domain of 
q m and by the IMS formula 

j=i Jr2 + 

The IMS error is easily seen to be bounded from below by some negative constant 
times .B~ 1 / 3 ~ 2<: ||V'|| 2 . By Proposition [ 



q«W > (e + CoB- 1/3 )|^ 2 |r, 

where we can make the constant Cq as large as we want, by choosing B large (using 
the properties of the support of X2,b)- 
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We turn to q m [^>i], and note that, with A^g from Appendix IA.21 

qmfoW > J r2 Ai, e (^Co + ^ 1/6 (Ci + y))^!! 2 + B- 1 ^\d p yj 1 \ 2 drdp. 

Since <; < 1/12 we can Taylor expand the first eigenvalue Ai,g to get a constant 
C > such that, on the support of ipi, 

Ai, e (co + 5- 1/6 (Ci + y' 

> Ai, e (Co) + ^ e (Co)(Ci + y) _ C (Ci + y)V 1/2 

> e + 5 (Ci + ^) V 1 ' 3 - C£T 

We insert this into q m to get 

dmW > e ||Viir+s- 1/3 / Is^l'+Jofci+^lViI'drd^cs-^+^iiVill 2 . 

Next, we use the modified Montgomery model from Appendix IA.3I to estimate the 
integral above, 

,2\2 



d, 



■M 2 + + y) foM 2 dTdp > A ljjCf (&)||ifc 



2 



We choose q = 1/48 and put the pieces together to obtain (|4.ip . 

We continue with the proof of (B). It is enough to prove (|4.2p for Q m (B), i.e., 

A2, Qm (B) > QoB + (70 + 7 )S 2/3 + 0(B 7/12 ), as £ -» oo, 

The inequality for Tt m (B) is then a direct consequence of Lemmas 12.51 and 12.81 
Let i/j^ and ip^ denote the first two normalized eigenfunctions of Q m (B) and 

assume that m is such that ^2,Q m (B) < ©o + wi?~ 1 / 3 for some w > 0, otherwise 

there is nothing to prove. 

We use the same cut-off function < Xi,b(p) < 1 as in the proof of (A), but 

with <; = 1/72. We also introduce the quadratic form q^ with the same action as 

q TO , but with an additional Dirichlet condition at \p\ = IB'' . For simplicity, we 

extend functions in the domain of q^ by zero for \p\ > 2_B S . We also denote by 

Qm(B) the corresponding self-adjoint operator. 
We start by showing that 

A 2 ,Q m ( B )>A 2!Q z, (B) +0(B- 00 ), asB^oc. (4.3) 
Let us write tpj — Xj-B^ k \ j, k — 1,2. By the IMS formula, it holds that 

2 

By Proposition ^. 11 we have for k = 1,2, as B — ► oo, 

||^||=0(B-~)||^||, q m [4 k) }=0(B-°°)\\^\\ 2 , and (4.4) 

B^il I |xUp)| 2 -|^ (fc) | 2 dTdp = 0(B— )ll^ (fe) H 2 . 
By the min-max principle we have 

i/'6span{i/'<' 1 ),j/' <2) } IIV^II 2 Q ^ 



A 2 , Qm(B) = max max q m [a^ (1) + /?^ (2) ]. (4.5) 



M" + |/3|"=l 
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Using Proposition 13.11 and (|4.4p we see that 

q™M (1) + /# (2) ] = qmM^ + Wf'l + 0(B-°°), as B oo. 
In the right-hand side we can write q^ instead of q m . It follows from (|4.4p that 

\\aipP + (3^ ] || 2 = 1 + 0(B-°°), as B -» oo. 
Using the min-max principle for q^ we have 

, . &] q m [^ 1) +^ 2) ] 
A 2 qb(b)= mm max < max j-. j-. — 

> dimV=2</, G V HVII 2 a,/9 \\ a ib\ X ) + M 2 ) ||2 

|q| 2 + |^| 2 = 1 11 ri 11 

= max q m [a^ 1} + /3ip[ 2) ] + 0{B-°°), as B -> oo. 

a,/3 
l"| 2 + l/9| 2 = l 



Combining this with (|4. 5|) we get 

Next, we show the existence of a positive constant 7 such that the inequality 

\2,Qg(B)>®0 + (%+j)B- 1/3 + O{B- 5 / 12 ), asB-^oc, (4.6) 

holds for all m € Z for which is bounded. Let e = B~ 1,/3 and ^£> € Dom(q^). 
We write q®bP] as 

CM=e[ |^| 2 + (r + Co) 2 |V'| 2 dTdp 

r f ~ D— 1/6 , n 2 \\ 2 D-l/3 

+ (1 - s) jf |3 r «| 2 + (r + Co + TT7 (Ci + y) J l*| 2 + TT7 M * dp 



I^X.( Jr ' / "( S + T))'l*l'^ 



For sufficiently large B we use the support of ip and the assumption |£i| < C\ to 
bound the last integral, uniformly in m, 

T^I( B " ,6 ( f '4)) 3 '*' 2d ^ 



as B — > 00. 



We get that q^M satisfies 

q£M>e / + (r + Co) 2 |^| 2 drdp 

Jr 2 . 

+ (l-e) jf Ai l0 ^ + T 3 7 (Ci + y)jM a + T 3 7 IM drdp 

+ 0(B- 2 ' 3+ ^)Uf, asB^oo, 

where Ai.g is the lowest eigenvalue of the de Gennes model, see Appendix lA.2l We 
use that ip has bounded support and estimate, using the Taylor expansion of \\ } g, 
as B — > 00, 

/ R2+ A - (Co + ^(C 1 + ^))l^drdp 

" L ( 90+ (T^^ 1/3 (Ci + y) 2 )H 2 drdp + 0(B-V 2+6s) ||^|| 2 
> ^ (©o + ^^^-Vs^ + ^^l^pdrdp + OCB-V^^ip 
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In the last inequality we also used that (1 — e)~ 2 — (1 — e) _1 = 0(B^ 1 ^ 3 ) together 
with Corollary 13. 21 Inserting in q^j we have, with the choice ? = 1/72, as B — > oo, 

&\>e I |d T ^| 2 + (T + Co) 2 H 2 dTdp+(l- e )e / \yj\ 2 drdp (4.7) 
+ £ 2 \d P ^\ 2 + £o(Ci + y)V| 2 drdp + 0(B-^ 12 )U\\ 2 

= e{g{Q ) ® i)m + (i - £)e ||vil 2 + y^(i ® MCi))M + o(s- 5/12 ). 

Here the operators Q and were introduced in Appendix [21 We note that the 
variables r and p are separated in the last expression, so if we denote by T the 
operator corresponding to the form on the right-hand side above then we have 

Ai,t = £A 1)5 (Co) + (1 - e)e + A li5J r(Ci)B- 1/3 + 0(£r 5/12 ) 
= 6 + A li? (Ci)B" 1/3 + 0(B- 5/12 ), as B — oo. 

Denote by 

70 =A 2 ,£?(Co)-Ai,g;(Co) and 7^ = inf (A 2 ^f(Ci) - A 1 ^f(Ci)) 

the spectral gaps for the de Gennes model and Montgomery model respectively. 
Since both \x,g((o) and A x rf(Ci) are simple eigenvalues, and Ci is varying in a 
compact set, it follows that both and 7^ are strictly positive. 
For the second eigenvalue of T we get 

A 2 ,t = Oo + (A m7 (Ci) + min(7 e , 7j? ))B- 1 / 3 + (D(B~ 5 / 12 ), as B -> 00. 

If we choose 7=5 min(7g, 7rf), we see that 

A 2 , T > e + (A^Ci) + 7)B" 1/3 + 0{B- b ' 12 ) 

> ©o + (70 + 7)^ 1/3 + O(B- 5 / 12 ), asB^oo, (4.8) 

where in the last inequality we use that A 1 jTf(Ci) > A. ^(C) = 7o- By (|4.7|) it 
follows that A^qdjb) > Aj,t for all j, so by (|4.8[) we get (14. 6p . 

The proof of (|4"T2j) for Q m (B) now follows by combining ((43)) and l|4.6p . □ 

5. Calculating good trial states 

5.1. Statement. We provide three different estimates of (the lowest point in) the 
spectrum of Tt m (B). Each is superior to the others in a specific parameter regime. 

When combined with Proposition 14. 1[ Theorem 15. II will give two-sided bounds 
on the ground state energy X\.n m (B)- 

Theorem 5.1. Let So > be the constant from Lemma \A . 2\ There exist constants 
Co, Ci? (2, Aj, j = 0, ... ,5, and polynomials \4(8) = SqS 2 + A4, As(5) with constant 
term As, and Ag(C3) quadratic with quadratic term <5oCfj such that with 

6= (m-^-CoVS-CiB 1/3 )5- lAi -C2, and 
Cs = m - I - Co VB - Cii? 1/3 - C 2 S 1/6 , 
f/ie following holds: 
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(i) For any K 2 / 6 > there exist K 2 / 6 and B 2 i 6 such that if \S\ < K 2 i 6 B 1 / 6 



and B > B 2 / e then 



dist 



3 1 

AjS^/ 6 + A 4 (<S)£r 4 / 6 , Spec(-H m (i?) 

3=0 



<K 2/6 (l + \6\ 2 )B- 5 / e . (5.1) 



(ii) For any K 1 / e > there exist Ki/§ and B 1 / e such that if \S\ < Ki/ e and 
B > Bi/q then 



dist 



3 1 
J2^B- J/6 + \i{5)B-^ + X 5 (S)B-^ e ,Spec^-n m (B: 

3=0 



< K 1/6 B~ 6 / 6 . 

(iii) For any Kq > there exist Kq and Bq such that if \^\ < Kq and B > Bq 
then 



dist 



5 

]T XjB-^ 6 + A 6 (C 3 )B- 6/6 , S P ec(-W m (i?) 

3=0 



< K B- J / 6 . (5.3) 



We emphasize that the constants Xj, j = 0, . . . , 5 and Co, Ci ancl C2 agree with 
the constants in Theorem ll.il In particular Ao, Ai and A2 are given in (|1.6|) . 

We will spend the rest of this section to prove this theorem. We start by proving 
part (iii), which is most detailed, and then go back and implement the necessary 
modifications for (i) and (ii). 



5.2. Proof of Theorem HUJ (iii). 



5.2.1. Outline. We will use a trial function %jj on the support of which the operators 
Qm(B) and TL m {B) agree (after a change of coordinates), see Section [5.2.131 below. 
This implies that it is enough to prove the theorem for Q m (B) instead of Tt m (B). 

We start by expanding the operator Q m {B) in powers ofB" 1 / 6 . Then we use 
the Grusin method [TH [53] to produce a trial state that agrees with (I5.3[) . 

5.2.2. Expansion of Q m (B). We will write 

Cs = m - I - (qV~B - CiB 1/3 - C 2 S 1/6 , 

and once the optimal values of Co, Ci an(1 C2 ar e determined we will write those as 
Co, Ci an d C2 respectively and insert them in the definition of C3- 
We expand the operator f) := Q m (B) in the form 

00 

f) ~ *)jB- j/e . (5.4) 
3=0 

The expansion (|5.4[) is to be understood as follows: For any function / in 5(M^) 
and for all N it holds that 

i)f = J2 B- j/ %f + 0(5- (Ar+1)/6 ) 
3=0 

in the sense of L 2 (ffi.^_). We recall the formula (|1.25p for Q m (B) and expand each 
term in a Taylor series, valid for small values of B^ 1 / 3, p and B~ x I 2 t. We will need 
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the first seven operators, 



fjo = -a? + (T + Cor, 

2 

[ )l = 2(Ci + ^)(r + Co), 



2 



with 



[)3-2C 3 (r + Co) + 2C 2 (Ci + y) +t) 3 , (5 ' 5) 

2 

f)4 = 2C 3 (Ci + y) + Cl + 6 4 , 
[)5 = 2C2C3+r(3r + 4Co)C2 + b5, 

f)6 = Cl + t(3t + 4Co)C3 + 2((Co + r)p 2 + 2rCi)C 3 + fj 8 . 

^3 = 29 T + r(r + 2Co)(r + Co) 

1)4 = CoV + |(6Ci + P 2 )r 2 + 4Co (Ci + y ) r 

f) 5 = -2rd 2 p + ^(4Co + r) + 2Cip 2 (r + Co) + K\r 



H„ = P 3 P + 2r9 T + ^p 6 + |V + Cl V + jr 4 + 4 Co r 3 + 3C 2 r 2 



Remark 5.2. The operators f)j, 3 < j < 6 are chosen to be independent of £2 and 
£3. For future reference, we also note that the linear terms in C2 in the operators f)3, 
f)4 and f)5 re-appear as linear terms in £3 in the operators t)4, f)5 and t)6 respectively. 

We study asymptotic expansions on the form 

00 00 
X^^XjB-^ 6 and ip(r, p) ~ ^ -(r, p)B~ i/6 . 
i=o i=o 

We want to find A., (as small as possible!) and ipj such that 

(f)-A)V'~0. (5.6) 

The expansion (|5.6p is to be understood term- wise, i.e., 

k 

52(t) i -\j)il>k-j = 0. (5.7) 

To start the Grusin approach, we need a function to project on. A study of (|5.7p 
for k — provides us with that. 

5.2.3. Order _B°/ 6 , a starting point. At this order equation (|5.7p reads 

fto - A )V'o = 0. 

Notice that f)o does not act in the p variable. We let tpo(T,p) — uq(t)(Pq(p). We do 
not want <po to be identically equal to zero, so we are led to solve 

-u% + (t + Co) 2 u q = A u , u{,(0) = 0. 

This is the eigenvalue equation for the well-known de Gennes operator £7 (Co) from 
Appendix IA.2I The smallest eigenvalue Ao is simple and given by 

A = e 

which is obtained for 

Co = Co = £0, (5.8) 
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see Lemma IA.2I The eigenfunction uq is positive and belongs to <S(]R+) (see 
Lemma [A. 3p . We may also assume that uq is normalized, J Q u^dr = 1. 

5.2.4. Higher orders, the Grusin approach. In this subsection we will implement 
the Grusin method [TH [23] which provides us with a systematic way of calculating 
a trial state for t) := Q m {B). We start by introducing some notation. First, we let 

oo 

= f) - (l)o - A ) ~ J2 ^ B ~ j/6 + A °' 

and notice that 

oo 

Srj-A-^fe-A^i^/ 6 . 
J=i 

Next, we introduce operators i?+ : L 2 (M) -> L 2 (K^), R~ : L 2 (M? + ) -> L 2 (R) and 
Eq : L 2 (R^_) — > L 2 (K 2 h ) (and in the corresponding Schwartz spaces) as 

(R + <P)(t,p) = <p(p)uo(T), 



(R-.f)(p) = / /(T,p)uo(T)dT, and 
Jo 



E Q =I® i?reg- 

Here i? reg is the regularized resolvent introduced in ljA.2|) . We introduce the matrix 
operators TL and £oj both acting in the Hilbcrt space L 2 (S? + ) L 2 (M), as 

n - *) , ^ & - (* « + 



By noting that 

%-Xq R a 
R- 



So 



we see that 



Notice that /C = 0(£? - 1 ' 6 ) as an operator on the Schwartz space. 
Let N <EN. Then by (5~ 



N 

HSoYsi-^^ = I + (-1) N JC N+1 = I + 0(B-( N +V/ 6 ), asS^oo. (5.10) 

3=0 

By (|5.9p we see that K) is given by 

K i = MM - X)E y [(61) - X)E } j -H6i) - X)R+^ 

Let us dehne E N , E N '+, E N ~ and E N >± via 



,E N >- E 

Then, by (j5"TU)) 



n .— Zn.± ):=£ ^- 1 ) 3Ki - 

' 3=0 



(t) - X)E N + R+E N - = 1 + 0(B- {N+1 V*) 
ft - X)E N '+ + R+E N - ± = 0{B-^ N+1 ^ 6 ), 
R-E N =0(B- ( - N+1 ^ 6 ), 

R- E N,+ = 1 + ( B -(iV+l)M 



(5.11) 
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as B — > oo. Assume that 

iV 

<p(j,) =Y J <Pj{p)B- j/6 , with E N ' ± ^ = 0(B- ( - N+1 ^ 6 ), asB^co, (5.12) 

3=0 

with ipj 6 5(R) for all j. Then by inserting the vector (0, <p) in (|5.10p and using 
the second formula of (|5.11[) we find that 

(\)-\)E N -+ V = 0{B- l - N+1 ^ 6 ), asB^oo. (5.13) 

We expand 

N 

E N <+ = J2 E t B ~ 3/& + 0(fl- (JV+1)/6 ), and 



3=0 
N 

E N '± = E EfB-V* + <D(B-( N+ V/ 6 ), as B -> oo, 

3=0 



(5.14) 



where the operators 



(E+f)(T,p)= E i-^(U E ^-kj)u (r)f(p) (5.15) 



ii,...,i»e{i,...,i} m=i 
/i+...+/i=j 



and 



(25f /)(/>) - / «o(r) 



3-1 

i=o h,...,i i+ ie{i,...,j} 

h+-+h+i=j 



x E (i)i - X h+1 ) 



U (T)f( P )dT (5.16) 



are independent of N. Here we use the convention that the summand is just \j) 
for i = 0. We write 

E N >+<p = E E+ Vk B-^' & + 0(i?-( A '+ 1 )/ 6 ), as B -> oo, (5.17) 

0<j,k<N 
j+k<N 

and define for j 6 {0, . . . , TV} the function to be the coefficient in front of B^^ 6 

in the sum on the right-hand side in (|5 . 1 7[) . We use (|5.15p to get a formula for ipj , 



^(r,p)= E (-!)*( II Mh m - Mj)uo<p k - (5.18) 

it!,...,/ 

We define our trial state by 



/ /.•• ! : j\ m=i 

h+...+k+k=j 



N 

i,{r,p)=YUr,P)B- j ' & . (5.19) 

3=0 

Since the operators involved are continuous (uniformly in B > 1), considered on 
Schwartz functions (the functions ijjj will be Schwartz functions, see Section l5.2.12p . 
it follows from (|5.13[) that 

(t]-\)iP = 0(B- (N+1)/6 ), asy3->oo. 

In particular l|5.7p holds for all < k < N. 
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Before we start with the calculations, let us note that the condition (|5.12p on 
the right reads that 

k 

Ef<p k -j = 0, for all < k < N. (5.20) 

3=1 

We also introduce the notation 

E t = E t~ X r 

5.2.5. Order B" 1 ^ , calculation of \\. To calculate Ai we need the operator E 1 . 
By (|A.3[) we see that, for any function / S L 2 (R), it holds that 

Eff = -RTfa - X 1 )R+f = -R- Ufa + y) (t + Co) - A X W/ = Ax/, 

and so E^ is just a multiplication operator. The equation (|5.20[) reads for k = 1 

= E 1 ip Q = \itpo- 
Since we do not want ipo to be identically zero, we find that 



Ai = and Ef = 



as an operator. 



5.2.6. Order B 1 / 3 , calculation 0/A2. Using that Ef — 0, the equation (|5.20[) for 
k = 2 reads 

= -E£<po = -R-[fomi - (f)2 - A 2 )]i?+^ 

2 2 

i(Ci + y ) (r + Co)i?rc g 2(Ci + y) (r + Co) 



OO 

"0 



(i 



5 2 p 2 \2 



</3 uo dr 

(-^ + (l-4. 1 )( Cl + ^) 2 -A 2 )^. 



(5.21) 



The number fci is defined in (|A.4[) and satisfies (1 — 4fei) = <5o by (|A.5[> . We choose 
A 2 = A 2 (Ci) to be the smallest eigenvalue of the operator —^2 + <5o(Ci + ■• and 
optimize in Ci , using the analysis in Appendix IA.31 to get 

A 2 = 2- 2 / 3 i> «5 ( ; /3 = 70, and Ci = Ci = C = (2<5 )" 1/3 C (5.22) 

From now on ipo denotes the corresponding eigenfunction, normalized in L 2 (R). 
We note that this also fixes the function ipo(T, p) = uo(t)¥>o(p)- 

5.2.7. Order B- k / 6 , k > 3. On the level B~ k l & , our unknowns are Afc, <fk-2 and 
also, up to some level, C2 and C3- The following procedure will determine tpk-2 and 
Afc. For tpk-2 to exist, it must be possible to solve 

k 

-E±ip k - 2 = ^Eftpk-j. 

3=3 

As we saw above — E 2 is the Montgomery operator minus its lowest eigenvalue. 
Thus, we want the solvability of the differential equation 

o2 2 9 \ ^ 

+<5o(Ci + yV-A 2 W2 = ]>>3 ± ^-3- ( 5 - 23 ) 



dp 



2 



3=3 
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But this is equivalent to the condition that the right hand side is orthogonal to the 
ground state tpo, i.e. 

J2E?Vk- j )dp = 0. (5.24) 

From (|5.24[) we get a formula for the unknown Afe , 

,k-\ 

\j=3 



/oo 1 \ 
^° E t^k- 3 +E±ip a J dp. 



Note that the right hand side will at some levels depend on £2 and £3. As soon as 
it does, we will minimize Afe over £2 and £3. When Afe is determined we find tpk~2 
by inverting the Montgomery operator in (|5.23|) . 

We calculate the first terms, the coefficients that appear are introduced in Ap- 
pendix El 

5.2.8. Order B^ 1 / 2 , calculation 0/A3. For A3 we get 
<p (E£<po) dp 

= {tpouo, [l) 3 - f)i£o(f)2 - A 2 ) - (t)2 - A 2 )£ fh + i)iE hiE i)i]u tpo) ■ 

Calculations, using the formulas (|5.5[) for fjj and the calculations in Appendix |A1 
give 

7 1 

((fiouo, f)3?Wo) = - g«o(0) 2 + Co - 2^0' 

(ipoua, t)iE (f)2 - A 2 )zwo) = 0, 
(ip u , (h 2 - \2)E a t) 1 u <p ) = 0, 
(ip u , i)iE \)iEoi)iu cpo) = 8fc 2 M 3 , 
which implies that 

A3 = -^o(O) 2 + - + 8fc 2 Af 3 . (5.25) 

5.2.9. Order B~ 2 / 3 , calculation of A4. For A4 we get 

/oo 
<Po{Eftpi + E^(f ) dp 
-OO 

= {tpouo, [({j 3 - A 3 ) - fji£ (f)2 - A 2 ) - (f) 2 - A 2 )£ f)i + fci-Efofh-EbfiijuoPi) 
+ (ipoUQ, [f) 4 - f)i-Eo(f)3 ~ A 3 ) - (f) 2 - A 2 ) J Bo(t) 2 - A 2 ) - (f) 3 - A 3 )£ f)i 
+ f)i^ f)i^o(f)2 - A 2 ) + f)l^o(t)2 - A 2 )S t)i + (t) 2 - A 2 )Eof)l£ f)l 
- f)i£'o^i-Bof)i-E , o{)i]uo^o)- 

Instead of calculating all these integrals explicitly we only look for their dependence 
on (2 and £3- It turns out that they are all independent of £3 and that £2 appears in 
the following integrals, which we calculate with help of the relations in Appendix IA1 
(we let const denote any constant independent of £2)' The quadratic dependence 
on £2 is given by 

(tp u 0l t) 4 ?Wo) = C2 + const, 
-(<poU , (rj 2 - A 2 )£'o(f)2 - A 2 )tWo) = -4fei(f. 
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(5.26) 



(5.27) 



The linear dependence on £2 is given by 

(ipou , (f) 3 - A 3 )u </?i) = 2C 2 A/oa + const, 
-((p u ,f)iE (f) 2 - \ 2 )u tpi) = -4M 1 1 fciC2, 
-(ip o u , (f) 2 - A 2 )£ f)itWi) = ~4M 1 1 fciC2, 
(tpouo, t)iE t)iE (b 2 ~ A 2 )tWo) = 8M^ fc 2 C2, 
((p u , t)iE (t) 2 - \ 2 )E i)iu ipo) = 8Mg i0 fc 2 C2 + const, 
((p u , (f) 2 - \2)E t) 1 E t) 1 u (p ) = 8M^ k 2 C 2 . 

The result is 

A 4 = (1 - 4/c 1 )C 2 + K M li( l - 4fc i) + 12M o,ofe)C2 + const 
= <5oCl + 2(M 1 )1 <5 + l2Ml Q k 2 )& + const. 
We see that A4 depends on £2 as a parabola and is minimal if 

C2 = C2 = ~ {Ml ,60 + l2M$ Q k 2 ) = -M X l - MloUk^ 1 . (5.28) 

This fixes the value of A4. 

5.2.10. Order B~ 5 / e , calculation 0/A5. For A5 we get 

/oo 
ip (Ef(f 2 +E±tpi+ Efipo) dp. 
-00 

To calculate all the integrals corresponding to A5 in full detail would be too cum- 
bersome for the presentation. The importance of this step is the dependence on 
(3. Thus we only calculate the integrals involving £3. First, we get some integrals 
involving £3 and £ 2 . 

(<^ "0, f)5lW0> = 2C2C3, 

-(ipaUQ, (h 2 - A 2 )-E (f)3 - X 3 )u Q ip ) = -4fciC 2 C3 + const, 
~(ip u (t) 3 ~ A 3 )-E (f)2 - A 2 )u Vo) = -4fciC2C3 + const. 



Here const denotes a constant that is independent of (3. By Remark 15.21 the linear 
terms in £3 (independent of £2) are the same as the linear terms in £ 2 for A4, 
see (|5.26[) . The result is that A5 is given by 

A 5 = [2A/ 1 1 (1 - Akx) + 2(1 - 4fci)C2 + 24M 2 /c 2 ]C 3 + const. 

By the choice of £2 in (|5.28p we see that the coefficient in front of £3 is zero, so A5 
is independent of C3 ■ We continue with the calculation of the dependence of £3 in 
A 6 . 

5.2.11. Order B^ 1 , calculation of \§. We do not calculate \§ in its full detail. The 
important part is its dependence on £3. There are two integrals that give rise to 
quadratic terms Q: 

(<PoUo, t)eUo(po) = C3 + lower order terms in £3, and, 

-((p u , (rj 3 - A 3 )£ , o(t)3 - X 3 )u ip Q ) = -4fciCf + lower order terms in C3- 

The result is that A$ depends on £3 quadratically with coefficient 1 — 4fci = 8q > 
in front of Cf- We introduce the constants £ 3 and C via 

A 6 = A 6 (C3) = £o(C3-C 3 ) 2 +C. (5-29) 
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5.2.12. Order B~ 7 / 6 and _B~ 8 / 6 , regularity properties. We note in order to obtain 
the functions c^ 5 and tp^ we should continue the calculations to the scales B^'^ 
and B -8 / 6 respectively. 

We end this section by noting that the facts that uq and (fo are Schwartz functions 
(Lemmas IA.3I and IA.8|) and that the resolvents of the de Gennes and Montgomery 
operators both maps the corresponding Schwartz space continuously to itself (Lem- 
mas IA.3I and IA.9|) . imply that the functions if>j, j — 0, . . . , 6, given by (|5.18[) all 



belong to S (Mi) 



- ' 



5.2.13. End of proof of Theorem \5.1\ (Hi). The calculations above provide us with 
Xj, and functions tjjj G S(R+), j — 0, ... ,6. We note that among the constants Aj, 
A6 is the only one that depends on (3, as in (|5.29p . Moreover, by carefully following 
(3 through the calculations of ipj, we find that ipo, 01 and ip2 do not depend on £3, 
while 

03 = 03,0 + 03, 1(3 

04 = 04,0 + 04,l(3 + "04,2(3 



05 = 05,0 + 05, 1(3 + "05,2(3 + 05,3(1 



06 = 06,0 + 06,lC3 + 06,2(1 + 06,3(3 + 06,4(1, 

where all the involved functions belong to <S(R+). We let \b be a usual smooth 
cut-off function, satisfying 

Xb(t, P ) = 1 on {(r,p) | < r < ^B 1 ' 2 , \p\ < f 5 1/3 }, and (5.30) 

supp( Xs (r,p)) c{(r,p) | < r < l -B 1 ' 2 , \p\ < \b^), 

and with Neumann condition at r = 0, (<9 t x_b)(0, p) = 0. With 

6 6 6 

A = 5^A J -B- J '/«, $( T ,p) =XBj2^ B ~ i/6 ' and Qm = J2 i >i B ~ j/6 + R 7 

j=0 j=0 3=0 

we can write 

6 k 6 6 

(^-a^^t^-a,^^^ 

k=0 j=0 k=lj=k 

+EEfe- A i)( 1 -^)^^ /6 +EEfe- A i)( 1 -^)^+6-;5-( w )/ 6 . 

fc=0 j=0 k=l j=k 

The first sum vanishes according to (|5.7p . Since |^3 1 < K the second sum and i?70 
are both bounded by a constant (independent of (3) times B~ 7 / 6 . The last two 
sums are of order 0(B~°°), since all functions 0j G <5>(R+). We therefore get the 
existence of constants Kq and Bq such that 

||(£L-A)0|| <K B- 7 / 6 \m 
for B > Bq. By the spectral theorem we conclude (|5.3I) . □ 



5.3. Proof of Theorem 15.11 (i). We use the same approach as in the proof of (iii). 
Actually, we repeat the same calculations, but with the differences that (3 = and 
(2 = C'2 + S. The result is that the operators f)o and f)i are independent of 5, while t)2 
involves S linearly in the form 2<J(t + (o), f)3 involves S linearly and f)4 quadratically. 
By keeping track of the S in the calculations, we find that A will have the form 

A = A + A 2 £- 2 / 6 + A 3 B- 3/6 + A 4 (<5)-B~ 4/6 , (5.31) 
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where Ao = 9o, A2 = 70, A3 is given in (|5.25[) and A4(<5) = A4 + S0S 2 where A4 is the 
coefficient calculated in Section 15.2.91 and So > is the constant from Lemma IA.2I 
The trial state has the form 

= XB (Vo + ipiB- 1/6 + (0 2 ,o + 2 ,i£)£T 2/6 + (0 3 ,o + 03,i<5 + ^2S 2 )B-V 6 

+ (0 4 , o + fa,i8 + 04, 2 <5 2 + 4 , 3 <$ 3 )5~ 4/6 ) , (5.32) 

where all involved functions belong to <S(R?_). 

We write Q m (B) = Ylj=o fyjB~^ 6 + R5 and organize the terms in {Q m (B) — X)ip 

as 

44 44 
(Qm(B) A)V7 = - X^-iB-"'* + EE^' - WwB-W* 

k=0 j=0 k=l j=k 

4 4 

+ R ^ + EEfe - w - xB^k-jB-*' 6 

k=a j=o 

4 4 

+ E E^ - A ^(! - x B )^ + 4-,i?- (fc+4)/6 . 

fe=l j=fc 

The first double sum is zero by (|5.7[) . For the second one, we use (|5.32[) and (|5.31[) 
to bound it in L -norm by (recall that \S\ < K 2 / 6 B^ 6 ) 

C"(l + \5\ 2 )B- 5 ^ 6 (l + \8\B- 1 ' 6 + \S\ 2 B- 2 / & + \S\ 3 B-' 3 / 6 ) < C"(l + \S\ 2 )B- 5 / 6 . 

Since the original operator depends quadratically on 5, it follows that 

||Jfe$| < (i + |^| 2 )s- 5/6 ||^||. 

Again, the last two sums are 0(B~ CO ), since the involved functions belong to the 
Schwartz space. We get, with A as in (|5.3f |> - that there exist constants K 2 /e and 
B2/6 such that 

||(Q m -A)0|| <£ 2/6 (l + |5| 2 )£r 5 / 6 ||0|| 
for B > i?2/6- By the spectral theorem we conclude (15. f p . □ 

5.4. Proof of Theorem 15.11 (ii). We repeat the Grusin calculation from the 
previous step once more, but this time we take one more term in the expansion. 
The result is a function of the form 

= XB (Vo + 0l£T 1/6 + (02,0 + Tp2,lS)B- 2 / 6 + (03,0 + 03,1<5 + 3 , 2 £ 2 )£r 3/6 
+ (0 4 ,o + 04,1-5 + 4 ,2<5 2 + 04, 3 5 3 )£T 4/6 
+ (05,0 + 05,1 5 + 5 ,2<5 2 + 05,3<5 3 + 05,4<5 4 )£r 5 / 6 ) 

where all involved functions belongs to 6>(R+). We also get 

A = A + A 2 B- 2 / 6 + A3B- 3 / 6 + \ A {S)B- A I & + X 5 (S)B- 5 / 6 , (5.33) 

with the same Ao, A 2 , A3 and A4(<5) as in (|5.3ip . Moreover, A5 (5) depends on S as 
a polynomial with constant coefficient A5, calculated in Section [5.2.101 

Xs(S) = A 5 + aiS + a 2 S 2 + a 3 S 3 

for some ai, a 2 and (13 in TSL. 
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We write Q m (B) = Ylj=o + ^ 6 an< ^ organize the terms in (Q m — X)ijj as 

5 5 5 5 

(Q m (B) A)^ = X)Efe - ^)^-i s ~* /6 + EE^ - A,)y, fe+5 _,i?-( fc+5 )/ 6 

k=0 j=0 k=l j=k 

5 5 

k=0 3=0 

+ E E(^ - A i)(! - Xs)^ + 5^i?- (fc+5)/6 . 

fc— 1 J — 

Again, the first double sum is zero. The second two terms are of order B~ 6 / 6 , 
uniformly for bounded 8, and the last two sums are of order 0(B~°°). We get, 
with A as in (|5.33|) . the existence of constants K\/q and B\/q such that 

||(Qm-A)^|| <^ 1/6 S- 6 / 6 ||^|| 

for B > Bx/g. We use the spectral theorem to conclude inequality (|5.2[) . □ 



6. Refined lower bounds 

6.1. Statement. We combine the lower bounds from Section [5] with part (i) and 
(ii) from Theorem l5.1l to have the improved lower bound. 

Theorem 6.1. Let K > 0. With the constants Co? Cij (,2, Aj, j — 0, . . . , 5 and 

8q, and with (3 from Theorem \ 5.1\ there exist constants Bq and Kq such that if 
IC3I ^ Ko an d B > Bo then 

Ai,w m ( B ) > A i? + X 2 B i/6 + \ 3 B 3/6 + A 4 S 2/6 + X 5 B 1/6 + K. (6.1) 

6.2. Proof. We divide the proof into several parts, depending on the size of (3. 
Actually, for most values of £3 we prove stronger results. We remind the reader 
that Ao = Oo and A2 = 70 ■ 

6.2.1. Proof for K 1/2 B 1 / 2 < |£ 3 |. It follows from Lemmas l2~5l and l2~7l that there 
exist constants if 1/2 and B X j 2 such that for all m £ Z, \^\ > K 1 / 2 B 1 / 2 , B > Bi/ 2 
and for all tp that satisfies Q m (B)ip = \{B)ip with X(B) < O B + ujB 2 / 3 it holds 
that 

q m [V>] > (e + i)IHI 2 , 

which clearly implies (|6.1[) . 

6.2.2. Proof for K 1/3 B^ 3 < | Cs I < K 1/2 B^ 2 . Assume that |f 3 | < ^i/ 2 ^ 1/2 and 
that B > B 1 / 2 . By Lemmas 12. 101 and !2.1 H it follows that there exist constants if 1/3 
and -B1/3 such that for all mgZ, |^| > K 1 / 3 B 1 ' 3 and B > B 1 / 3 it holds that 

q m m > (eo + ^o + l)^- 1 / 3 )!^!! 2 , W>eDom(q m ). 

By possibly changing if 1/3 and -B1/3 slightly, it follows from Lemmas 12.51 and 12.81 
that for all m € Z such that \^\ > K1/3B 1 ' 3 , B > B X j 3 and for all tp that satisfies 
Q m (B)ip = X(B)ip with X(B) < Q B + ujB 2 ' 3 it holds that 

q m m > (e + ( 7 o + l)B- 1/3 )IH| 2 , 
from which (|6.1[) follows. 
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6.2.3. Proof for K 5/16 B 5 / 16 < |£ 3 | < K ll3 B x l z . Assume that | Cs I < K 1/3 B^ 3 and 
that B > Bx/3 with the constants if 1/3 and Bi/ 3 from the previous step. We define 
6 by the equation 

m=| + CoS 1/2 + (Ci+^ 1/3 , 

where Co and Ci are the constants from (|5.8[) and (|5.22[) . and note that the condition 
ICs I < K 1/3 B^ 3 implies that \5\ < C for some constant C. From Lemma IA.7I it 
follows that for all |<5| < C there exist a positive constant C pos such that 

\m(Ci +S)> X lM (Ci) + C pos \S\ 2 . (6.2) 

For \S\ > CiB^ 1 / 48 we combine (|6T2]l with Proposition O (A) to find that 

q m [V] > (e +7oS- 1/3 + (C pos C\ 2 - C)B^ 8 )U\\ 2 . 

for sufficiently large B. If we choose Ci sufficiently large, we get the existence of a 
positive constant C' pos such that for C\B~ X I 48 < \S\ < C and for all B large enough 
it holds that 

q m l4>} > (Qo+joB-^ + C^B- 3 /*)^. 

Finally, we invoke Lemmas 12.51 and 12.81 decrease the constant C' pos slightly to 
Cp OS if necessary, to get existence of positive constants K 5 /i§ and -B5/16 such that 
if K 5/16 B 5 / 16 < |Ca| < K 1/3 B 1 ^ 3 and B > B 5/16 , then for all ip that satisfies 
Q m (B)il> = X{B)ip with X(B) < Q B + ujB 2 / 3 it holds that 

q™M > {Q B + j B 2 / 3 + C; os B 5 ^)U\\ 2 . 

This inequality is also stronger than (|6.1I) . 

6.2.4. Proof for K^qB 1 / 6 < | Cs I < K 5/16 B^ W . Assume that |C 3 | < K 5/16 B 5 / 16 
and -B > -B5/16 with the constants K^ng and -B5/16 from the previous step. We 
introduce S as 

m = § + Co5 1/2 + CiS 1/3 + (Ca + <^ 1/6 , 

where Co , Ci and £2 are the constants from (|5.8p , (|5.22p and (|5.28fl respectively, and 
note that \5\ < CB 7 / 48 for some constant C. 

Since 7/48 < 1/6 we may apply Theorem 15.1 f i). It follows by Proposition l4.il 
that 

^h,n m (B) = A + X 2 B- 2 / 6 + A3B- 3 / 6 + X 4 (S)B-^ 6 + 0{{l + \S\ 2 )B- 5 / 6 ), 

as B — » 00. 

For large \S\ and £?, we have 

X 4 (S)B- 4 ^ + 0((1 + \5\ 2 )B- 5 ' 6 ) > A 4 + 2 

where A4 is the constant in (|5.2T[) . We let K ± / 6 and -B1/6 correspond to the constants 
for which Kx/qB 1 / 6 < \( 3 \ < K 5/w B r °/ 16 and B > B 1/6 implies that 

Xi, Hm (B) > XqB + X 2 B 4 / & + A 3 £ 3/6 + (A 4 + 1)£ 2/6 . 
This clearly implies (|6.1[) . 
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6.2.5. Proof for K < \( 3 \ < K^qB 1 ^ 6 , final step. Assume that |Cs| < Ki/ e B 1/6 
and that B > £>i/6, with the constants if 1/6 and -B1/6 from the previous step. 
Again, we let S be given by 

m = § + CoS 1/2 + CiS 1/3 + (C2 + S)B^. 

This time \5\ < C, for some constant C, so we can apply Theorem I5.1f ii) com- 
bined with Proposition 14.11 Recall that \4(S) = 5q5 2 + A4 and Ag(<5) = A5 + 
a\5 + 02<5 2 + a3(5 3 . We rewrite the last two terms in the eigenvalue expansion from 
Theorem 15. lf ii) as 

A 4 (<5)B- 4 / 6 + \ 5 (S)B- 5 / 6 

= A4B- 4 / 6 + A5B- 5 / 6 + 5 6 2 B- 4 / 6 + ( ai 6 + a 2 S 2 + a 3 S 3 )B-^ 6 . 

If we choose C sufficiently large and C > CB^ 1 / 6 then the term 5o8 2 B~ i / e will 
dominate both (<zi<5 + a 2 S 2 + a 3 5 3 )B~ 5 / 6 and the error which is bounded by some 
constant times £>~ 6 / 6 , indeed, we can get that all these three terms are bounded 
from below by \5 C 2 B- Q / Q . 

Therefore we find that there exist constants Bo and Kq such that if it holds that 
K < |Cs 1 < K 1/6 B 1 / 6 and B > B a then 

Ai,« ra (B) > A S + A 2 B 4 / 6 + A3B 3 / 6 + A4S 2 / 6 + A5B 1 / 6 + K. 
This finishes the proof of Theorem 16.11 □ 

7. Proof of Theorem 11.11 

For bounded (3 we combine Theorem I5.1f iii) with Proposition 14.11 to get the 
asymptotic formula 

Kn m {B) = ^oB + A 2 B 4 / 6 + A3B 3 / 6 + A 4 5 2/6 + A 5 B 1/6 + A 6 (C 3 ) + 0(£T 1/6 ) (7.1) 

as B — > 00. Comparing the lower bound from (|6.ip with (|7.ip . we find that the 
lowest eigenvalue is smallest for bounded £3, and that its asymptotic expansion 
then is given by (|7.ip . For bounded C3 we see from (|5.29[) that the smallest value 
of Ae(Cs) is given for £3 = (3. However, since m must be an integer, we are not free 
to choose ( 3 = (3(111, B) arbitrarily. With 

A B = inf \( 3 (m,B)-( 3 \ 

as in (|1.4[) we find that the smallest possible Ae(Cs) is given by 

\ 6 = S Q A 2 B +C. (7.2) 
This finishes the proof of Theorem 11.11 □ 



8. MONOTONICITY OF Xi,n(B), PROOF OF THEOREM 11.31 

We first note that by perturbation theory it holds that 

\,-H(B),+ - K,h(b),- 

for all B > 0. 

From Theorem 1 1.1 1 we know that the lowest eigenvalue X\.u(B) °f satisfies 
Ai,h(b) ^eQB + X 2 B 2 / 3 + \ 3 B 1 / 2 + \ 4 B 1 ^ + X 5 B 1 ^ + 5oA 2 B +C + 0(B- 1 / 6 ) (8.2) 
where 

A B = inf |C 3 (m,S)-C3| 
and Cs(m, B) = m - § - ( Vb - (1B 1 / 3 - (2B 1 / 6 . 
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It is proved in [5] that the derivatives A^ r yt B \ + satisfies 



liminf X' ln(B) + ^ hmsup - liminf (A 1 . H(s+e) - A 1/H(B) ) 

B — >-oo v 1 e >o+ £ " — >OQ 

limsupA^ n(B) < liminf -lira sup (A 1H(B) - X in r B _ e) ) 



We start with the right derivative X' ± 

X-l,h{b+ s ) ~ Ai,w( B ) g(B + e) 
— Uo H 



and use 



to write 



9(B) A| +e -A| , f(B 
r oo - 



(8.3) 



/(£) 



e £ e 

where 

= A 2 B 2 / 3 + A3S 1 / 2 + A 4 B 1/3 + AgB 1 / 6 + C 

and f(B) is a function satisfying lims^oo /(B) = 0. For any fixed e > we clearly 
have 

lim g(B + s)-g(B) = 0. 



A 2 —A 

Consider the term involving B+e — - 
rriB+s such that 

A B = 

A B+e = 
We note that 



We note that there exist integers mj and 



m B - - - ( Vb - (iB 1/3 - C2B 1 



/6 



C3 



and 



B 



m B +e 



a/6 



A B = inf 



< 



m - I - Co VB - CiB 1/3 - C2B 1 / 6 - ( 3 
m B+£ - ^ - (o\^B - CiB 1/3 - C2B 1 / 6 - Cs 



Using this, and the fact that < A B < 1/2 for all B, we get by the triangle 
inequality 



A 2 _ a 2 

^B+e B 



"(A 



B+E 



> 



1 ?N /B+7-VS 7 {B + e) 1 ' 3 - B 1 ^ 7 (B + e) 1 ' 6 - B 1 / 6 

7T + Co h Cl 1" C2 

2 £ £ £ 



The right-hand side tends to —1/2 as B — ► 00, so for any fixed £ > we get 

iliminf(A 2 B+E -A|)>-i 

£ S— >oo z 

Inserting these calculations in (|8.3[) it follows that 



liminf X' lH , B) + > 6 - -S . 

According to (|1.7p the right-hand side Qq — ^Sq > 0. This finishes the proof of (|1.9p . 
The same calculations give (|1.10l) for the left-derivative A^ 



H(B)- 



. We conclude the 



proof of Theorem ll.3l bv noting that the equations (|1.9p and (|8.ip imply that Ai^rm 
is increasing for large B. □ 

Appendix A. Model operators 

In this appendix, we consider two self-adjoint model operators. The first one 
is an operator in L 2 (K + ) that was introduced by Saint- James and de Gennes [2"2"] . 
The second one is an operator in L 2 (M.) first studied by Montgomery in [19] . 
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A.l. A general Lemma. We start by giving a general lemma that will enable us 
to give moment formulas for the two operators under study. 

Lemma A.l. Let — oo < a < j3 < oo, and p G C 1 [a, f3\ . (If a = — oo and (3 = +oo 
then we assume that \im x ^±r X1 p{x) = +oo). Assume that for some A € 1 and 
u G L 2 (a,/3) it holds that 

—u" + pu = Xu for all x G [a, 0\. 

Then, for any polynomial b, it holds that 

[b'" + A(X-p)b' - 2p'b]u 2 dx = [2b{uf + b"u 2 - 2b' W + 2(A -p)6u 2 ] (A.l) 

Proof. In the case a — — oo and/or — +oo, the additional assumption on p 
implies that u decays exponentially at a and/or (3 (the proof is the same as in 
Lemma IA.8|) . 

One could use the same reasoning as in pQ. However, a simple calculation shows 
that the derivative of the expressions inside the brackets in the right-hand side 
of (|A.1[) equals the integrand on the left-hand side. □ 

A.2. The de Gennes operator. For £ G K we define the operator G(£) as the 
self-adjoint Neumann extension in L 2 (M. + ), acting as 

(G(0U)(X) =-u"(x) + (x + OMx): 

u'(0) = 0. 

Denote by \i,g(£) the lowest eigenvalue of G(0 and 6o = inf^gR Ai i g(£). We refer 
to [5] for a discussion of the results summarized in the following lemma. 

Lemma A.2. The function £ i— ► \i,g(£) is smooth. Moreover, 

(1) lim^ +co Ai,£f(£) = +oo. 

(2) Bm € __ oo Ai, (0 = 1. 

(3) The function Ai i g(£) attains its minimum value Oo, \ < ®o < 1, at a 
unique point £o < 0. 

(4) Ai : £f(£) is decreasing for £ < £o Gf&d increasing for £ > £q. 

(5) XTie number Sq :— ^A" g (£o) satisfies < <$o < 1. 

If we denote by i>£ the normalized eigenfunction of corresponding to the 
eigenvalue \i,g(£), then we introduce the regularized resolvent 

(A.2, 

(0 .9 || U£, 

and let i? rcg = -R rC g(£o) an d u o = W| . 

Lemma A. 3 ( 6 , Lemma A. 5). The function uq belongs to <S(R+) and R re g maps 
<S(R+) continuously into itself. 

Lemma A. 4 ([\\, equations (2.34)-(2.36)). The following equalities hold 

oo />oo 

u^dx = 1, / (x + £,o)uq dx = 0, 

( A -3) 

( x + t ) 2 u 2 dx='f, J {x + Co) 3 uldx=^f>. 
We introduce the integrals 

p CO 

kj(0= (x + 0u[R ieg (0(x + 0Yudx (A.4) 
Jo 

and kj = fcj(£o)- 
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Lemma A. 5 ([BJ, Proposition A. 3). It holds that 

So = lXi,gfa) = 1 - 4*!. (A.5) 

Remark A.6. Numerical calculations of the constants £o, ©o, #o and uq(0) has 
been carried out in [2J. 

We give a new approach. It is readily seen that the decaying normalized solution 
to Q(£)u = Ai,g(0« is given by u{x) = ce-^ x+ ^ 2 H i (Ai (x + f). Hcrc 
c denotes a normalization constant and iij, is the Hermite function, that solves 
—y"(x) + 2xy'(x) — 2vy(x) = 0. The boundary condition lt'(0) = transforms into 

(Ai, e (0 - 1 ) H i(\i, a {Q-3)it)-€ H ttx 1 . {t)-i)($) = and since Q o = Kafa) = Co, 

we find that £o should be the largest (it is negative!) number that solves 

(£ 2 -l)%£.-3)(0-£%e»-l)(0=0. 

Numerical calculations in Mathematica give 

Co « -0.76818365314, 6 « 0.59010612495, 
tt o (0) w 0.87304313851, 5 « 0.58551290029. 

A. 3. The Montgomery operator. Next, we turn to the Montgomery operator 
M.(Q, C S K, defined as the self-adjoint operator in L 2 (R) acting as 

(M(C)«)(P) = + (C + P 2 ) 2 "W, -oo < p < oo 

Denote by Xi,m(C) the lowest eigenvalue of M.(Q and Vq = inf^gR Ai^ (£). 

Lemma A. 7 f |19L I21L I13j). The function C, h- > Xi.m(C) * s smooth and satisfies 
limf —>±oo Ai,x (C) = +oo. Moreover, the minimal value vq > o/ Ai^CC) * s 
attained at a unique point £ < 0, and A'/ jW (C) > 0. 

Let us denote by <£> the eigenfunction corresponding to Auvi(C)- It is known that 
such an eigenfunction belongs to C°° (R) . We show that ip and its derivatives decay 
exponentially, which implies that ip belongs to <S(K). 

Lemma A.8. Let ip be the ground state of M(£). For any < a < 1/3 and 

nonnegative integer k there exist a constant Ck such that 

J e 2a ^ 3 (M 2 + |pV| 2 + l(5p) fe ^| 2 ) dp < C k J M 2 dp. (A.6) 

Proof. Let a < 1/3 be given. For e > we define v s (p) — (|p|/(l +e|p|)) 3 . Then, 
for fixed p, v £ (p) is monotonically increasing to |p| 3 as e — * 0. Moreover it holds 
that 

K(p)\ < 3|p| 2 (A.7) 
for all e > and all pel. We let Xx,M and X2,m denote the same functions as in 
the proof of Lemma 12.121 The IMS formula gives 

/ \d P ( X 2,MZ av ^)\ 2 + f(C + P 2 ) 2 - Ai,m(C)) \X2. M e av ^\ 2 dp 
Jr v 7 

= Jjd p ( X 2,M^M 2 dp (A.8) 
We use the Cauchy-Schwarz inequality on the right-hand side, to get, for any ? > 0, 

/ \d p ( X 2,M^M 2 dp 
Jm. 

< (1 + f \av' E e av ^X2,Mv\ 2 dp +{! + -) [ | {d pX 2,M)e av ^\ 2 dp (A.9) 
Jr ? Jr 
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We recall that \3 p X2,m\ < I2/M for all p and that d P X2,M has support in the set 
{p E K. I M < p < 2M}. We implement (fA~7j) and (fO|) in fX8|) to find that 

\d P (x2, M e aV ^)\ 2 + ((C + P 2 ) 2 - Ai,m(C) - (1 + 09a 2 |p| 4 ) \x%M^^ dp 

<(i + ^ e a(2M)3 / R M 2 dp (A.10) 

Since a < 1/3 we can choose ? so small that (1 + s)9a 2 < 1. With this choice of <r 
we can find M so large that 



(C + p 2 ) 2 -A liA4 (C)-(l + 09a 2 H 4 ) 



> 1 



for all p on the support of X2.M- This together with the trivial bound for small p 
settles the result for ip. We might also use the first term in (jA.lOjl to prove the 
result for d p . The statement for p k ip follows from this by decreasing a (or to be 
more precise, prove the result for Lp for 1/3 > a' > a and then decrease this a' to 
a). The result for higher derivatives is now a consequence of induction, using the 
eigenvalue equation. □ 

Let us define the regularized resolvent i? reg as 

~ / (M(o - A^torv u ±<p, 

RregU = < (A.ll) 

10 u || ip. 

We show that if u and its derivatives decay exponentially, then the same is true 
for R Teg u. 

Lemma A. 9. Let < a < 1/3. Assume that u, p k u and (d p ) k u belong to 
L 2 (R, e 211 ^ dp) for all non-negative integers k. Then, for any b < a, R reg u, p l R reg u 
and (d l p )R reg u belong to L 2 (R, e 2b ' p ' dp) for all non-negative integers I. 

Proof. Let w = i? rcg u, so that 

- d 2 p w + (C + P 2 ) 2 w - \ ltM (()w = u, (A.12) 
with u as in the assumptions. Let e > and let v e be the same function as in the 
proof of the previous lemma, v £ (p) — (\p\/(l + e\p\)) 3 . We also let xi.m and X2,m 
be the same cut-off functions as in the proof of Lemma 12.121 An integration by 
parts gives 

I d p ( X 2,Me^ w)\ 2 + \{( + p 2 ) X2, M z aVc w\ 2 dp 



= Km{Q / \X2,M^ aVe M^P+ / \d p ( X 2,M^)w\ d p+ / \ X2tM e av *\wudp. 

We use the Cauchy-Schwarz inequality for the last term and move terms to the 
left-hand side, 

|5 p (x2,Me QU ^)| 2 +P(p)|e^u;| 2 dp<2 / \ X 2.Me av *u\ 2 dp, 

Jr 



where 



P(p) = (|(C + P 2 )x2,m| 2 - (Ai,a4(C)-2)| X2 ,m| 2 - \d P { X 2,M)+av' M i\ 2 ). 

By (|A.7[) , the first term in P is dominant for large p, so if we choose M large enough 
we have P(p) > 1 for all p on the support of X2,m, and thus we get, for such M, 
that 



'{|p|>Af} 



|<9p(x2,j\/e a ' Ue w)| 2 + \e av 'w\ 2 dp < 2 / | X 2,Me Q1, ^| 2 dp. (A.13) 

JR 
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The right-hand side is clearly bounded by 2 JL |u| 2 e 2a ' p ' 3 dp which is bounded by 
assumption. We let e — > and use monotone convergence to conclude 

/ \w\ 2 e 2a \p\ 3 dp < 2 f \u\ 2 e 2a ^ 3 dp. 

J{\p\>M} JR 

The estimate for \p\ < M is trivial. This proves the statement in our lemma for 
w (with b = a). The estimate for p l w is simple if we just decrease a to b, so 
that p l is dominated by the exponential exp(2(6 — a)|/o| 3 ) for large p. For the 
first derivative d p w we might use (|A.13[) and for higher derivatives we continue by 
induction, using (|A.12jl . □ 

We will several times encounter a variant of the Montgomery operator. For 
k > 0, we denote by A4(() the operator 

(M(Qu)(p) = -u"(p)+k(c+^) 2 u(p), -oo<p<co (A.14) 
A change of coordinates 

transforms into 

M(0 = k l ' 3 2- 2 ' 3 M{k l / 3 2 l / 3 Q 
and so for the lowest eigenvalue it holds that 

\m(0 = fcV3 2 -2/3 Al)M(A; l/3 2 l/3 C) _ (A.15) 

In the case when k = 6q we write 

C = (26 )- 1/3 C 
and we get that A 1 r^(C) 1S minimal for C = C an( i that 

X lM (0 = 2~ 2 / 3 6 1 /3 \ 1 MO = 2- 2/3 ^ /3 ^o = %. (A.16) 
For the case k = 6q we also introduce the moments 

f°° t~ o 2 \ l 

M lk = J (C + yJ <Pj<Pkdp. 

Here ipo is the first normalized eigenfunction of for k = 5q, and tpj, j > 1 are 

constructed via the Grusin method in Section [5l 

Lemma A. 10. It holds that 

< = 1, M 1 ,o = 0, Ml = —^—, and M$ „ - ' 



3(2<5 ) 2 / 3 ' °'° 6S 3(2<5 ) 2 / 3 ' 

Proof. The first formula is just the normalization of <po and the second one follows 
by a perturbation argument, just as for the de Gennes model. 

We use Lemma [A. II to calculate M l for I > 2. Indeed, with p(p) = (C + 
and A = 70 , the formula (|A.1|) becomes 



^dp = 0. (A.17) 

The choice b{p) = p 21 ^ 1 gives M^ 1 provided that the previous moment formulas 
are known. Especially the choices b(p) = p and b(p) = p 3 give the announced 
formulas for Mq and Mq . □ 
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